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Abstract 

A Bose gas in an external potential is studied by means of the local density 
approximation. Analytical results are derived for the thermodynamic proper- 
ties of an ideal Bose gas in a generic power-law trapping potential, and their 
dependence on the mutual interaction of atoms in the case of a non-ideal Bose 
gas. 
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One of the most striking consequences of quantum statistics is the condensation of an 
ideal Bose-Einstein gas where the zero momentum state can become macroscopically oc- 
cupied at a sufficiently low temperature [[[]. Some of the features of the transition to a 
superfluid phase exhibited by the Bose fluid 4 He is often interpreted essentially as a result of 
this Bose-Einstein condensation (BEC), the strongest degeneracy effect of a Boson system 
||. For many years, it was considered hopeless to experimentally observe BEC in an atomic 
gas with weak interactions. With the development of techniques to trap and cool atoms, 
BEC was recently observed directly in dilute atomic vapors |3|-[|]. The new experimental 
achievements have stimulated great interest in the theoretical study of inhomogeneous Bose 
gases. 

The thermodynamic properties of trapped atomic Bose gases undergoing BEC can be 
altered by the spatially varying trapping potential. The interaction between atoms may have 
a significant effect on the thermodynamic properties. There have been several investigations 
analyzing the dependence of the critical temperature on the trapping potential and weak 
interaction in the Bose gas The thermodynamic properties of Bose gases in an external 

potential have also been discussed in Ref. |]], but no analytic results have been given. Here 
we shall derive some analytical expressions for the thermodynamic properties of a trapped 
non-interacting and interacting Bose gas under the local density approximation. 

Let us first consider the case of an ideal Bose gas trapped in an external potential. When 
the energy level spacing due to the trapping potential is much smaller than the thermal 
energy kT = /3 _1 , the local density approximation is adequate p|-[TT[]. Space may then be 
divided into small cells, and in each cell we may consider the trapping potential V(r) to be 
constant. With a trivial modification to indicate r dependence explicitly we may directly 
extend the formula of the grand potential, which can be found in any textbook on statistical 
mechanics, e.g. Ref. to write the local grand potential above the critical temperature T c 
as 
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where p is the chemical potential, the thermal wave length A = (2irh 2 /mkT^j 2 , and 
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From this grand potential the local density of the Bose gas in the external potential is 
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This is the same as [p.l|| . For the generic power-law potential discussed in Ref. || 
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where Li, Li and L3 are the linear sizes of the volume, and E\, E2, £3, p, I and q the 
parameters, integrating the local grand potential ([I]) over the whole volume, we find the 
total grand potential of the system to be 

n BE (T,p) = -kT\- 3 J g 5/2 (()dr 
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where z = e^, rj — 1/p + l/l + 1/q + 1/2, and 
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In the derivation we have used the definition of the function g u (x) and the formula for the 
gamma function 
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From the total grand thermodynamic potential the average atom number can be obtained 



as 
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This is in agreement with the result in Ref. || obtained directly from the local density p(r) 
instead of the grand partition function. 
The total entropy can be obtained as 
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Although, for obvious reasons, the volume and pressure are not useful thermodynamic vari- 
ables for a trapped gas ||, we may still consider an analogue to the capacity for the ho- 
mogeneous system at constant volume with a fixed number of atoms, and define this heat 
capacity as 

c ( THr(§) N , (io) 

which includes the work done against the potential as the energy of the gas is changed. 

So far we have discussed the case for a temperature above T c . When the temperature T 
is below the BEC critical temperature T c , the chemical potential /x = 0, and only the normal 
component of the Bose gas has a contribution to the heat capacity. From Eqs. (§) and ( [TDD 
we can then obtain the heat capacity (T < T c ) to be: 
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As a limit, at the critical temperature T c , the above expression reduces to 

C(T-) = {r] + 1) ^ + fl 9 ^ 1) (N)k. (12) 

When the temperature T is above the critical temperature T c , at a fixed iV the chemical 
potential jj, is a function of temperature T. In order to obtain the heat capacity, we have to 
calculate the derivative of \x with respect to temperature T while holding iV fixed. Using 
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from Eq. (§) we find 
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From Eqs. (|T0| ) and (|TJ) the heat capacity above T c is then given by 
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At T = T c + , the chemical potential /i = 0, so the above expression reduces to 
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By comparing with Eq. fll2l) , the discontinuity of heat capacity at the critical point is 



(17) 



AC = C{T~ 



C{T+) = iv + 1)2 . 9 ^ 1) (N)k. 



(18) 



From the Eqs. (jTTD , (]17|) and fljjf ), we see that the heat capacity depends only on the particle 
number and the exponents p, I and q of the external potential. We can further examine the 
continuity of C(T) and its derivatives. As we know, g v (l) is divergent for 77 < 1. So, the heat 
capacity becomes discontinuous at the critical temperature T c when 77 > 1. From Eqs. (|i~T|) 
and (|16"1), the derivative dC(T)/dT is discontinuous at T c for either 77 < 1 or 77 > 1. This 
is in agreement with Ref. 0. From the analytical expressions we also see that the heat 
capacity at the critical point grows with increasing 77. 

For a rigid wall box, 77 = 1/2. It is straightforward to check that the above expressions 
indeed reduce to the known results in textbooks of statistical mechanics (e.g. Ref. (]]]). As 
a special case let us consider the harmonic potential 
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We find the heat capacity C(T C + ), C(T C ) and AC to be 
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This agrees with Refs. 

Along similar lines we may analyze the thermodynamic properties of a trapped interact- 
ing Bose gas. In the mean field approximation, the single particle effective potential can be 
written as: 



V eff (v) = V(r) + 2a\ 2 p(r), 
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where a is the s-wave scattering length, which characterizes the strength of interaction 
between two atoms. We assume \a\ C A for the weak interaction considered here. By 
similarity with Eq. ([]]), for the trapped interacting Bose gas the local grand potential can 
be written as 
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where 
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Expanding the function 05/2(0 in Eq. ( f24|) according to the small parameter a/ A, and keeping 
terms only up to the lowest order in a/A, we have 
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which coincides with the first order result of the virial expansion of the grand partition 



function [13], This is not surprising since the mean field approximation is accurate to the 
first order. From the grand potential we can obtain the local density 
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which, compared with the density of an ideal Bose gas, means the effective fugacity is 
exp[— (3(V(r) + 4aA 2 p(r) — //)]. However, from a naive consideration based on the effective 
interaction one would think that it should be exp[— /3(V(r) + 2aA 2 p(r) — //)] instead. The 
factor two in the atomic density expression naturally appears by taking the derivative of the 
grand potential with respect to \i. 

As the counterpart to Eq. @, we can find the total grand potential for the interacting 
gas at T > T c to be 
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where 
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When a = 0, Eq. fl28|) recovers expression (||) for the non-interacting Bose gas, as it should. 
The average number of atoms is now given by 
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where we have used 
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As we know, g v+ i(z) is a bounded, positive, monotonically increasing function of z, and 
Fs,v,r){%) is positive. From Eq. (|30| ) we see that the trapped interacting Bose gas has a higher 
chemical potential than an ideal Bose gas when the atomic number and temperature are the 
same in both systems. 



From Eq. (|30D for (N 1 ) we can calculate the derivative dfi/dT at a fixed iV 7 , and the 
result is 
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Corresponding to Eq. (Q), the total entropy of the trapped interacting Bose gas can then be 
derived as 
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Calculating the derivative of the entropy (S^), we obtain the heat capacity above the BEC 
critical temperature to be 
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The last two formulas show the correction due to the mutual interaction of atoms to ther- 
modynamic quantities such as the entropy and heat capacity. 

The above derivations for the trapped interacting Bose gas is adequate only for the gas 
phase (without BEC), and furnish no information whatsoever about the condensed phase, 



or the nature of the condensation [IB]. Any attempts to use the same equations in the 



regime of BEC will lead to some mistakes. According to Ref. []T3|] the mean field effective 
potential should be amended to include the contributions of atoms to Bose condensation. 



If we choose an appropriate effective mean field potential similar to Eq. (28) of Ref. [13 



we can, in principle, extend the above discussion to calculate thermodynamic quantities 
below the critical temperature. However, the derivation will be so involved that no simple 
analytical formulas can be expected to exist any more. 
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